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En ; ABSTRACT 

We have obtained Hubble Space Telescope/STIS low-resolution ultraviolet spectra of the X- 
ray pulsar 4U 1626-67 (=KZ TrA); 4U 1626-67 is unusual even among X-ray pulsars due to its 
ultra-short binary period (P=41.4 min) and remarkably low mass-function (< 1.3 x 10~ 6 Mq). 
The far-UV spectrum was exposed for a total of 32ks and has sufficient signal-to-noise to reveal 
numerous broad emission and prominent narrower absorption lines. Most of the absorption lines 
are consistent in strength with a purely interstellar origin. However, there is evidence that 
both C I and C IV require additional absorbing gas local to the system. In emission, the usual 
prominent lines of N V and He II are absent, whilst both O IV and O V are relatively strong. 
We further identify a rarely seen feature at ~1660A as the O III] multiplet. Our ultraviolet 
spectra therefore provide independent support for the recent suggestion that the mass donor is 
the chemically fractionated core of either a C-O-Ne or O-Ne-Mg white dwarf; this was put forward 
to explain the results of Chandra high-resolution X-ray spectroscopy. The velocity profiles of the 
ultraviolet lines are in all cases broad and/or flat-topped, or perhaps even double-peaked for the 
highest ionization cases of O; in either case the ultraviolet line profiles are in broad agreement 
with the Doppler pairs found in the X-ray spectra. Both the X-ray and far-UV lines are plausibly 
formed in (or in an corona just above) a Keplerian accretion disc; the combination of ultraviolet 
and X-ray spectral data may provide a rich data set for follow-on detailed models of the disk 
dynamics and ionization structure in this highly unusual low-mass X-ray pulsar system. 

Subject headings: accretion, accretion disks — line: profiles — binaries: close — pulsars: individual 
(4U 1626-67, KZ TrA) — ultraviolet: stars — X-rays: stars 

1. INTRODUCTION 



1 Based on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science The X-ray Source 4U 1626-67 is One of the rare 

Institute, which is operated by the Association of Universi- cases i n w hi C h an X-ray pulsar is a member of a 

NAS f 5 r 2655r Ch " Astr ° n ° my * InC '* Undcr NASA COntraCt low mass X-ray binary (LMXB) system. Its 7.7s 

pulsations were first detected in 1977 by SAS-3 
(Rappaport et al. 1977), but even with further 
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extensive observations the pulse timing has yet 
to reveal any indication of a binary orbit, plac- 
ing very tight constraints on the mass function 
(< 1.3 x KT 6 M Q , Levine et al. 1988). However, 
following the identification of a faint, blue opti- 
cal counterpart (KZ TrA) by McClintock et al. 
(1977), extensive fast optical photoelectric pho- 
tometry provided convincing evidence for bina- 
rity. In addition to the direct optical pulses due 
to the reprocessing by the disc of the pulsed X- 
ray flux, Middleditch et al. (1981) uncovered the 
Dopplcr-shifted signal from optical pulses due to 
reprocessing on the donor star's heated face. This 
detection has now been confirmed (with the same 
instrumentation) by Chakrabarty (1998), and re- 
cently with a fast frame-transfer CCD camera 
(Chakrabarty et al. 2001). 

The derived orbital period of 41.4 min places 
4U 1626-67 as one of the ultra-short period 
LMXBs, i.e. those with P or h < 80 min, below 
the limit for a low-mass main sequence hydrogen- 
burning donor. In consequence there has been 
much controversy regarding the nature of the sec- 
ondary and the evolutionary history of the sys- 
tem. The earliest proposed donor types were a 
low-mass (0.08 M Q ) severely hydrogen-depleted 
and partially degenerate star (Nelson et al. 1986), 
or an even lower mass (^0.02 M ) helium or C-0 
white dwarf (Verbunt et al. 1990). Later moder- 
ate resolution X-ray spectra from ASCA (Angelini 
et al. 1995) and BeppoSAX (Owens et al. 1997) 
revealed emission line structures identified as due 
to Ne/O, indicating an overabundance of these el- 
ements in the system. To explain the excess Ne (a 
by-product of He burning) Angelini et al. (1995) 
proposed that the donor could be a a low-mass he- 
lium burning star, under-filling its Roche lobe and 
transferring material via a powerful stellar wind. 
The most recent progress has come from high res- 
olution X-ray grating spectroscopy from Chandra. 
Schulz et al. (2001) (hereafter SCH01) report on 
further abundance anomalies required to explain 
the strength of various absorption edges. From 
the inferred local abundance ratios they now ar- 
gue that the mass donor is most likely the 0.02 M 
chemically fractionated core of a C-O-Ne or O-Ne- 
Mg white dwarf. An evolutionary scenario has also 
been developed by Yungelson et al. (2002), indi- 
cating a probable 0.3-0.4 M Q C-0 white dwarf as 
progenitor. 



The Chandra spectra not only resolve the indi- 
vidual X-ray emission lines, but they also indicate 
double-peaked emission, especially in the case of 
Ne X A12. SCH01 interpret these as Doppler pairs 
arising most probably in or near the Keplerian disc 
flow. Indeed, this would be inconsistent with a 
more massive 0.08 M Q donor, for which i < 8°; 
whereas for a 0.02 M© white dwarf the maximum 
allowed i = 30° makes it much more plausible. 

4U 1626-67 has now been studied extensively in 
the X-ray regime and also photometrically in the 
optical. However, little useful spectral information 
has been extracted in the optical, as LMXBs gen- 
erally show only a few weak lines there. Cowley 
et al. (1988) obtained two spectra of 4U 1626-67 
which show only weak C III/N III 4650/60 emis- 
sion. In the far- ultraviolet (FUV), however, the 
spectra of LMXBs arc normally characterized by 
numerous strong emission lines due to the ioniz- 
ing flux of the incident X-rays. Hence, we have 
utilized the excellent UV capability of the Space 
Telescope Imaging Spectrograph (STIS, Woodgate 
et al. 1998) on-board the Hubble Space Telescope 
(HST) to obtain time-resolved spectra of 4U 1626- 
67 in both this line-rich FUV region and also in 
the near-UV (NUV). In this paper, we report our 
results on the time-averaged spectra, concentrat- 
ing on the FUV lines, and discuss our findings in 
the light of the recent X-ray spectral results in 
particular. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

We obtained HST /STIS UV observations of 
4U 1626-67 on 1999 April 2, May 31 and June 1 
totaling 12 HST orbits; we also extracted a fur- 
ther 4 orbits of archival data taken a year ear- 
lier on 1998 April 23 as part of another program. 
All data were taken with either the low resolution 
G140L (FUV) or G230L (NUV) gratings, which 
give useful wavelength ranges of ~1160-1700Aand 
~1700-3150A, at 1.2A and 3.2A resolutions re- 
spectively. The FUV and NUV MAMA detectors 
were operated in Time- Tag mode, which is then 
effectively photon-counting 2 . All these data were 
reduced by the On-The-Fly Reprocessing system 
at STScI, which utilizes the best currently avail- 

2 Time-resolved broad-band UV lightcurves are presented in 
Chakrabarty et al. (2001) 
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Fig. 1.— Far-UV spectrum of 4U 1626-67 (=KZ TrA) from a total exposure of 31.2 ks with HST/STIS. No 
correction for the small reddening of this object has been made, but to aid clarity 3 pixel boxcar smoothing 
has been applied. Numerous emission and absorption lines are present, the most significant of which have 
been labelled. Emission lines are labelled from beneath, whilst absorption or coinciding emission/ absorption 
lines are labelled from above. Note in particular the broad and complicated emission profiles; these are 
discussed further in the text. 



able calibration files, software and data parame- 
ters. The products provided include 2-D spectral 
images created by accumulating the photon event 
list over the exposure during a given HST orbit, 
as well as the wavelength and flux-calibrated 1-D 
spectra extracted from these. We then also ap- 
plied additional wavelength corrections using the 
STSDAS WAVECAL routine (if not previously 
completed), which makes use of the wavelength 
calibration (arc) spectra taken with each science 
exposure. Lastly, we derived a fully time-averaged 



spectra in each band from all the data (weight- 
ing each HST orbit's data according to exposure 
time), which amounts to total exposure times of 
31.2 ks (8.7 hrs) in the FUV (and 10.8 ks in the 
NUV). 

3. DATA ANALYSIS AND RESULTS 

In figure 1 we present the complete FUV spec- 
trum (not corrected for the small reddening to- 
wards this object), with 3 pixel boxcar smooth- 
ing applied to aid clarity (and with la error-bars 
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combined accordingly). The long exposure pro- 
vides good signal-to-noise, and numerous promi- 
nent broad emission lines and narrower absorp- 
tion features stand out. By contrast, the NUV 
spectrum shows no very strong emission features 
and only absorption lines (from Mg II, Mn II, Fe II 
and Zn II) of strength consistent with interstellar 
origin; hence we will only make occasional further 
reference to results from this NUV spectrum in the 
remainder of this paper. 

Apart from the very strong absorption feature 
at ~1215A due to Lya, which has prominent 
damping wings, the results of Gaussian fits (using 
the full-resolution data) to all the prominent FUV 
lines (measured line centers, FWHM and equiva- 
lent widths) are presented in table 1. Almost all of 
the emission lines exhibit complex structures and 
in most cases emission and absorption lines coin- 
cide. We will discuss the fits to these more com- 
plex profiles in turn, in the following sections. For 
the remaining absorption lines we fit each line with 
a single Gaussian, with center, width and normal- 
ization as free parameters. In the case of blends, 
if the separation of the two (or more) lines was 
greater than the spectral resolution, we included 
a Gaussian profile for each. 

Similarly for each line contributing to the emis- 
sion/absorption line complexes we used uncon- 
strained Gaussians. However, for the components 
of multiplets additional constraints were applied. 
We required the ratio of the components' centers 
(effectively the spacing) to agree with that of the 
vacuum wavelengths, and the FWHM of each com- 
ponent to be the same. The absorption compo- 
nents are likely dominated by interstellar absorp- 
tion. Hence we left the ratio of the fluxes (e.g. 
the doublet ratio) free to allow for the probable 
range of optical depth for the absorbing clouds. 
However, attempts at leaving the ratios free for 
the emission too, indicated that the data arc in- 
sufficient to meaningfully constrain these values, 
and often actually produced unphysical results. 
Hence, in most cases we simply assumed the opti- 
cally thick limit for the emission lines and used the 
appropriate relative intensities as tabulated in the 
National Institute of Standards and Technology 
(NIST) listings 3 . Most observational determina- 

3 available at http://physics.nist.gov/cgi- 

bin/AtData/lincs_form 



tions, as well as the results of modeling, (see e.g. 
Ko & Kallman 1994), indicate that these resonant 
emission lines are typically optically thick. Nev- 
ertheless, for the Si IV and C IV doublets we also 
examined the optically thin limit (where the dou- 
blet ratio is 2:1 for the lower to higher wavelength 
components), but found no substantive differences 
(see the relevant section for details). We there- 
fore tabulate (and show figures) based on the opti- 
cally thick assumption. Lastly, we also tried twin- 
Gaussians to represent two separate blue- and red- 
shifted emission regions, i.e. as was done for the 
double-peaked X-ray lines by SCH01. 



4 



Table 1 

Parameters of Gaussian fits to FUV absorption and emission lines in 4U 1626-67 



Compo 
-nent a 


ID 


■^■vac 


AA 

(A) 


V 

(kms- 1 ) 


FWHM 

(kms- 1 ) 


Flux (10~ 16 
erg cm -2 s -1 ) 


EW 

(A) 


+ 


cm 


1175. 64 b 


-0.33 ±0.52 


-85 ± 135 


1440 ± 330 


10.6 ±3.4 


0.68 ±0.22 






J- J- oo . u u 


+ fl 25 


—30 ± 65 


QOfl _|_ 1 cn 


—3.9 ± 1.6 


26 ± 1 1 


_ 


Si II 


1190.42 


0.25 ±0.11 


65 ±30 


360 ± 70 


-6^5 ± L4 


0.45 ±0.10 


— 


Si II C 


1193.29 


-0.01 ±0.06 


-5 ±15 


470 ± 40 


-9.2 ± 1.0 


0.63 ±0.07 


— 


+S III C 


1194.06 


-0.78 ±0.06 


-195 ±15 








- 


NI b 


1200.13 


0.16 ±0.06 


40 ± 15 


610 ± 40 


-12.3 ± 1.0 


0.88 ±0.07 





SII 


1250.58 


-0.98 ±0.40 


-235 ±95 


850 ± 230 


-10.5 ±2.8 


0.76 ±0.20 


- 


SII 


1253.81 


-0.02 ±0.11 


-5 ±25 


460 ± 90 


-5.9 ± 1.3 


0.42 ± 0.09 




Si TT C 


1 260 42 


_0 98 + n 05 

U.ZjO _l_ \J.\JO 


_ftz _i_ i n 

\JO -J— J-VJ 


Kf.() _|_ on 












u.uz zc u.uo 


lOU ZIZ ±u 








- 


CI 


1277.41 


0.33 ±0.21 


75 ±50 


320 ± 110 


-1.9 ±0.8 


0.13 ±0.06 


- 


CI 


1280.33 


-0.54 ±0.80 


-125 ± 190 


380 ± 280 


-1.2 ±1.0 


0.09 ±0.07 




O T 


1 QHO 1 7 

louz.l / 


U.U4 zt U. 14 


in ± on 
1U zt oU 


t^nn _l 1 1 n 
DUU zt 11U 


1 9 A _1_ O 7 
— 1Z.U ± Z. 1 


n qi i n in 

U.ol zt u.ly 




q; tt 

Ol II 


1 Qfl/1 Q7 
10U4.O I 


n 1/1 j_n 11 

U.14 zt U. 1 1 


oU zt ZO 


QOH -4- 7ft 

oyu zt /u 


— D.o ± i.y 


n /ift _i_ n i o 
U.40 zt U.lo 


— 


C I 


1328.83 


0.44 ±0.23 


100 ± 50 


530 ± 170 


-2.9 ± 1.0 


0.20 ±0.07 




C II 


1334.53 


-0.01 ± 0.16 


± 35 


500 ± 50 


-12.0 ± 1.9 


0.85 ± 14 


— 


CII 


1335.71 


-0.01 ±0.16 


0±35 


500 ± 50 


-4.6 ±2.4 


0.32 ±0.17 


+ 


IV 


1338.61 


0.37 ±0.47 


85 ± 105 


2780 ± 380 


9.7 ± 1.4 


0.69 ±0.10 


+ 


IV 


1342.99 


0.37 ± 0.47 


85 ± 105 


4~1 r~ o t \ 1 o O (~v 

2780 ± 380 


6.5 ± 1.0 


0.43 ± 0.07 


4. 


IV 


1343.51 


n 37 ± 47 


as + 1 05 


1 OVJ _l_ fj OVJ 


9.1 ± 1.4 


66 ±010 




/"I T 

C I 


1328.83 


0.30 ± 0.18 


70 ± 40 


440 ± 100 


—2.2 ± 0.7 


n 1 p* i r\ nr 

0.15 ± 0.05 




C II 


1334.53 


-0.21 ± 0.20 


-45 ± 45 


400 ± 70 


-9.1 ± 1.9 


0.64 ±0.13 




C II 


1335.71 


-0.21 ±0.20 


-45 ± 45 


400 ± 70 


-6.6 ±3.0 


0.47 ±0.21 


+ b* 


IV 


1338.61 


-3.61 ±0.38 


-805 ±85 


770 ± 230 


2.5 ±0.9 


0.18 ±0.07 


+ b* 


IV 


1342.99 


-3.61 ±0.38 


-805 ± 85 


770 ± 230 


1.6 ±0.6 


0.11 ±0.04 


+ b* 


IV 


1343.51 


-3.61 ±0.38 


-805 ± 85 


770 ± 230 


2.9 ±1.1 


0.21 ±0.08 


+r* 


IV 


1338.61 


2.88 ±0.77 


645 ±170 


1850 ± 430 


6.2 ± 1.6 


0.44 ±0.11 


i ^* 
+r 


\J IV 


1 Q A o on 
lo4z.yy 


z.oo zt U. it 


ft a rr 1 i 7n 
040 zt 1 / U 


1 QKC\ J- A on 

looU zt 4oU 


4.U zt l.U 


n. 07 j- n c\i 
U.z / zt U.U ( 


-i-r* 

-\-r 


O TV 


1O4O.01 


9 qq _i_ n 77 

Z.OO I U.l 1 


040 zn ± i u 


1 QKC\ _!_ A on 

±oOU zn 4oU 


( . 1 Z!Z l .o 


n c;o _i_ n 1 ^ 
U.OZ ZIZ u.io 


+ 


V 


1371.29 


0.27 ±0.31 


60 ±70 


2820 ± 190 


49.1 ±4.3 


3.66 ±0.32 


+ 


U V 


1 Q71 on 


/i /in j- n oo 
— 4.4y zt U.ZZ 


AQA _|_ Cfl 

— yoU zt oU 


you zt loU 


14. Z zt Z.O 


1 f\A j. n in 
1.U4 zt U.ly 


+r* 


V 


1371.29 


2.75 ±0.30 


600 ± 65 


1680 ± 150 


28.9 ±3.0 


2.12 ±0.22 


+ 


V 


1371.29 


0.20 ±0.18 


45 ±40 


1690 ± 70 


491 ± 20 


36.1 ±1.5 




V 


1371.29 


-0.21 ±0.19 


-45 ± 40 


1610 ±80 


-447 ±21 


32.9 ± 1.6 


+ 


Si IV 


1393.76 


(A fixed) 




850 ± 460 


12.5 ±12.2 


0.93 ±0.91 




Si IV 


1393.76 


(A fixed) 




590 ± 130 


-18.4 ±8.0 


1.36 ±0.59 


+ 


OIV] 


1401.35 


(A fixed) 




3060 ± 180 


42.5 ±4.5 


3.17 ±0.34 


+ 


Si IV 


1402.77 


(A fixed) 




840 ± 450 


10.1 ±9.9 


0.76 ±0.74 




Si IV 


1402.77 


(A fixed) 




590 ± 130 


-12.0 ±5.2 


0.90 ±0.39 


+ 


Si IV 


1393.76 


2.79 ±0.37 


600 ± 80 


2000 ± 160 


24.3 ±2.8 


1.81 ±0.21 




Si IV 


1393.76 


-0.02 ±0.08 


-5 ±15 


520 ± 50 


-10.3 ± 1.2 


0.76 ±0.09 


+ 


IV] 


1401.35 


1.52 ±0.27 


325 ± 60 


680 ± 10 


13.4 ±8.4 


1.00 ±0.63 


+ 


Si IV 


1402.77 


2.81 ±0.37 


600 ± 80 


1990 ± 160 


19.4 ±3.7 


1.46 ±0.28 




Si IV 


1402.77 


-0.02 ±0.08 


-5 ±15 


520 ± 50 


-13.8 ± 1.3 


1.03 ±0.10 
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3.1. Fits to individual Emission/ Absorption 
line complexes 

3.1.1. C 111X1176 

This profile appears to be fiat topped, or even 
double-peaked. However, given the lower signal- 
to-noise at this short wavelength end of the spec- 
tral range, such details are well within the un- 
certainties and a simple single Gaussian fit is as 
appropriate as any more complex model. The re- 
sulting line center is consistent within the uncer- 
tainties with the weighted average of the vacuum 
wavelengths of the closely spaced multiplet com- 
ponents. 



3.1.2. IV triplet A1342 

In addition to the three emission components 
of the O IV triplet, one must also include absorp- 
tion lines of C I A1329 and the C II A1335, A1336 
doublet in order to fit the observed profile. Even 
with the 1000 km s _1 span of the three components 
of the multiplet, the breadth of this emission fea- 
ture requires a very large FWHM of 2800 km s -1 
for each separate component (see fig. 2a). A 
better representation may well be a two Gaus- 
sian double-peaked profile for each component, 
as was used for the SCH01 X-ray line spectra. 
This more complex model does provide a better fit 
(fig. 2b), with a broad (FWHM=800km S - 1 ) blue- 
shifted component (V = — 800kms _1 ), but an 
even broader (FWHM=1900kms _1 ) red-shifted 
component (V — 600 km s -1 ). 

3.1.3. O V X1371 

This singlet line, unaffected by any expected 
strong interstellar absorption lines should be the 
most clear cut case. However, the single Gaussian 
fit again yields a very large FWHM of 2800 km s _1 , 
as well as being clearly a very poor fit (see fig. 3a). 
Once again motivated by the double-peaked X-ray 
lines, we have attempted such a fit, which works 
well (see fig. 3b). This yields velocity shifts of 
-950kms _1 and 650kms~ 1 for the blue- and red- 
shifted components respectively, with reduced, but 
like O IV A1342, differing FWHM of 950 and 1700 
kms -1 . For a disk representation the FWHM 
should be the same; hence this implies that there 



may be additional red-shifted emission contribu- 
tions. A fit using one blue-shifted and two red- 
shifted Gaussian components is also possible with 
a similar significance. 



Our attempt to fit the profile with a single 
emission peak and absorption to account for the 
dip at ~1369A failed, although it is still possi- 
ble that in addition to multiple emission peaks 
there is some absorption. For instance, phase- 
resolved STIS FUV spectroscopy of Her X-l (Vr- 
tilek et al. 2001; Boroson et al. 2000) also revealed 
unusual features in the O V line profile. In ad- 
dition to the broad and narrow emission compo- 
nents, absorption was required for the fits in Her 
X-l, which moved in phase with the binary and 
therefore was interpreted as a local effect. We 
note that published STIS FUV spectra (phased- 
averaged) of the X-ray transient XTE J1859+226 
(Haswell et al. 2002) also shows a similar O V 
structure to 4U 1626-67, but these authors do not 
discuss the profile. 

3.1.4. SiIVX139l XI 403 doublet + OlVjXUOl 
blend 

This is the most complex of the emission/absorption 
profiles in our spectrum. Our most basic model 
is to assume broad emission components for each 
of the Si IV doublet lines and the O IV line and 
narrower interstellar Si IV absorption, with the 
doublet ratio (DR) constrained to the optically 
thick case (=1.25). We first tried the model fixing 
all the central wavelengths to their vacuum values, 
and then relaxed this constraint. In this instance, 
significant differences in the line parameters re- 
sulted. Although, in both cases the model can be 
made to fit the profile quite well (apart from fail- 
ing to delineate the drop in flux between the two 
Si IV lines at ~ 1398 A), physical consideration of 
the resulting parameters shows them to be ques- 
tionable. The fits favor large line fluxes in both 
the emission and absorption components, which 
then balance out to provide the profile required. 
In the fully constrained fit almost all the emission 
flux comes from a very strong (EW=680kms -1 ) 
and very broad (FWHM^SOOOkms^ 1 ) O IV] line 
(see fig. 4a), whilst in the less constrained fit a 
large 600 kms -1 velocity shift of the Si IV dou- 
blet is found, together with very high fluxes and 
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Table 1 — Continued 



Compo 
-nent a 


ID 




AA 

\ A ) 


V 

( km s ) 


1 IT T TTTA IT 

FWHM 

^ Km S J 


Flux (10 ib 

— 2 — l \ 

erg cm s j 


EW 
( A ) 


+* 


Si IV 


1393.76 


-0.55 d 


-120 d 


870 d 


I3.2 d 


0.97 d 




Si IV 


1393.76 


-0.23 


-50 


540 


-13. 9 


1.04 


+ b* 


IV 


1401.35 


-3.77 


-805 


770 


8.8 


0.65 


+ r* 


IV 


1401.35 


3.00 


645 


1850 


21.6 


1.63 


+* 


Si IV 


1402.77 


-0.55 


-120 


870 


10.5 


0.79 


_* 


Si IV 


1402.77 


-0.23 


-50 


540 


-9.5 


0.71 




Si II 


1526.71 


0.11 ±0.07 


20 ± 14 


380 ± 30 


-7.2 ±0.8 


0.56 ±0.06 


+ 


CIV 


1548.20 


0.26 ±0.15 


50 ±30 


1490 ± 120 


27.5 ±5.2 


2.20 ±0.42 




CIV 


1548.20 


-0.02 ±0.07 


-5 ±15 


430 ± 40 


-14.2 ± 1.9 


1.14±0.15 


+ 


CIV 


1550.77 


0.26 ±0.15 


50 ±30 


1490 ± 120 


24.8 ±4.7 


1.98 ±0.37 




CIV 


1550.77 


-0.02 ±0.07 


-5 ±15 


430 ± 40 


-9.9 ± 1.8 


0.79 ±0.14 




CI 


1560.91 


-0.46 ±0.17 


-90 ±35 


320 ± 90 


-3.1 ± 1.0 


0.25 ±0.08 




Fe II 


1608.45 


0.00 ±0.13 


0±25 


330 ± 80 


-5.9 ± 1.6 


0.51 ±0.13 




Cl b 


1657.22 


0.40 ±0.20 


70 ±35 


260 ± 90 


-2.5 ± 1.1 


0.23 ±0.10 


+ 


III 


1661.20- 


0.60 ±0.30 


110 ±55 


580 ± 110 


17.1 ±4.6 


1.63 ±0.44 




(7 comps) -1669.30 














Al II 


1670.79 


-0.24 ±0.07 


-45 ± 15 


280 ± 30 


-9.5 ± 1.5 


0.91 ±0.15 



Note. — The entries in the table are grouped (separated by horizontal spaces) according to which lines 
were blended together and therefore required simultaneous fitting. 

a Coding to clarify the various components, especially in the case of the more complex absorp- 
tion/emission structures: + indicates emission component; - indicates absorption component; when two 
emission components (double-peaked case) are present r and b indicate the red- and blue-shifted lines. 
Also, where more than one model has been fit to a profile, the * indicates the results of the preferred 
model. 

b Weightcd average of vacuum wavelengths for closely spaced (AA < 1.2A resolution) multiplet compo- 
nents quoted, and single Gaussian used for fit 

c Single Gaussian used to fit blend of two lines with AA < 1.2A spectral resolution. 

d No errors are quoted for the parameters for this fit, as they only represent a local minimum in x 2 
space, that for which the Si IV emission/absorption fluxes are smallest (and most plausible). 

c O III] has 7 multiplet components in the range 1661. 2A to 1669.3A as given by NIST, we used a 
Gaussian for each. The wavelength shift and FWHM apply to each component (constrained to be the 
same), whilst we quote the summed flux and EW only. 
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Fig. 2. — Gaussian line profile fits to O IV triplet A1342 emission, with absorption lines of C I A1329 (solid 
line) and the C II A1335, A1336 doublet (dashed), a) Single emission component for each of the O IV triplet 
lines, note the very broad profile required (FWHM=2800kms~ 1 ). b) Assuming blue- (solid) and red-shifted 
(dotted) components for each triplet line provides a marginally better fit, but with significantly broader 
red-shifted than blue-shifted lines. 



2000 km s -1 line widths (see fig. 4b). Using a 
DR of 2 for the optically thin limit had little 
effect apart from the absorption intensities ad- 
justing to a similar ratio to balance the emission 
once again. Clearly much more complex emission 
components are needed. Since there is evidence 
that O V and possibly O IV profiles are better fit 
by an asymmetric twin-Gaussian model, we have 
tested whether such an O IV profile can fit the 
semi-forbidden line here. We have constrained 
the components to have the same velocity shifts, 
widths and relative normalizations as O IV A 1342, 
and find that this does work reasonably well (see 
fig. 4c). We note, however, that a similar double- 
peaked profile would be too broad to fit the su- 



perimposed Si IV features as they are relatively 
narrow. 



3.1.5. C IV doublet A1549 

As the baseline model for this doublet struc- 
ture we assume the presence of the two doublet 
emission lines (with DR=1.11 for optically thick 
limit) and corresponding doublet absorption lines, 
together with absorption at 156lA, which we at- 
tribute to C I. As can been seen in figure 5, 
this model does provide a reasonable fit to the 
data, with only a small 50kms _1 shift (compa- 
rable to the uncertainty) in the emission doublet 
wavelength required. The strong C IV emission 
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1360 1370 1380 
wavelength (A) 
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wavelength (A) 



1390 



Fig. 3. — Gaussian line profile fits to O V A1371. a) With a single emission component a very broad line 
is required, but this is a noticeably poor fit due to the dip at ^1369A, b) double-peaked emission, with 
red-shifted component (dotted line) much broader than the blue (solid), providing a much better fit. 



is a common feature of LMXB FUV spectra, but 
the absorption is unusually pronounced. Check- 
ing the optically thin limit for the emission (with 
DR=2), we found that the absorption strengths 
are in fact insensitive to this ratio. Even try- 
ing double-peaked emission for each component of 
emission doublet cannot reduce the depth of the 
absorption lines required, since the shape of the 
line profiles constrains each emission component to 
be almost as broad as the single Gaussians shown. 
We will discuss the strong C IV absorption and 
the probable C I absorption further in §§ 4.1.1. 
Also, we note that the width/ velocity spread of 
the emission lines is only about half that of the 
O V and O IV. 



3.1.6. O III] multiplet A1661 - 1669 

This profile can be adequately fit with a sin- 
gle emission Gaussian for each of the multiplet 
components, together with Al II A1671 and yet 
more C I (A1657) absorption. But, as for the 
O IV A1342 and O V A1371 lines, a very large 
FWHM^2800kms _1 is indicated for each compo- 
nent by the fit, and the feature is overall unusually 
strong for an LMXB, as compared to, for instance, 
the Si and C doublets. 

4. DISCUSSION 

4.1. Absorption and emission strengths 
and the local elemental abundances 

Henceforth, in the cases of the complex emis- 
sion and absorption profiles, we only consider the 
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Fig. 4.— Gaussian line profile fits to the complex Si IV A1394, A1403 doublet (solid lines) + O IV] A1401 
(dashed) blend. Both emission and absorption components are included for the Si IV doublet lines, a) 
Constrained fit with line centers given by vacuum wavelengths. This results in broad and strong O IV] 
emission, with narrow Si IV emission, balancing only slightly narrower absorption lines, b) With the line 
centers now unconstrained (apart from the separation of the doublet components), the fit yields much broader 
Si IV emission redshifted by as much as 600 km s -1 , and narrow O IV] shifted to balance the Si IV A1403 
absorption in an unlikely manner, c) Motivated by the fits to O IV A1371 and O V A1342, the single O IV] 
line has been replaced with a twin-Gaussian profile (dashed and dot-dashed), with separation, FWHM and 
relative normalizations constrained to that found for O IV A1342. Essentially the results are similar to 
a), though the Si IV emission is now slightly blue-shifted relative to its absorption. None of these fits is 
physically very satisfying, and in all likelihood the true underlying emission profiles are even more complex. 



fluxes/EWs given by the physically reasonable fits. 
These particular fits are indicated by asterisks in 
table 1 to aid cross-reference. 

4-1. 1- Interstellar versus local absorption 

At first glance, the equivalent widths (EWs) of 
the many absorption lines (if interstellar in origin) 



would appear rather large given the low column to- 
wards this source, E(B-V) = 0.11±0.1 (SCH01). 
However, consideration of the results from the 
HST Quasar Absorption Line Key Project sug- 
gests otherwise for most of the lines. Savage 
et al. (2000) report on the Galactic interstellar 
features they identified in the spectra of this large 
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Fig. 5. — Gaussian line profile fit to C IV dou- 
blet A1549, including both emission and absorp- 
tion, together with probable C I A1561 absorp- 
tion. Fairly broad (FWHM=1500kms^ 1 ) and 
strong C IV is required, with much narrower 
(FWHM=430kms~ 1 ) absorption. 



sample of 83 quasars. Their sample includes a 
range of E(B - V) from 0.014 to 0.155, as deter- 
mined from radio H I measurements; therefore the 
range encompasses the reddening of 4U 1626-67. 
In table 2 we summarize our measured equivalent 
widths for the numerous absorption lines (includ- 
ing those superimposed on top of emission), and 
present the comparative results from Savage et al. 
(2000), both their tabulated median and rms val- 
ues for the EWs and an estimate of the maximum 
they measured (from their fig. 5). In almost all 
cases where a comparison is possible, our mea- 
sured values for 4U 1626-67 lie comfortably within 
the range for the quasars. The only exceptions are 
the Si IV A1393/1402, C II A1335/1336 and C IV 
A1548/1550. As noted above the Si IV A1394/1403 
lines are blended with O IV] making fitting to 
this profile rather uncertain. Given the fact that 



the other Si lines (all Si II) have EWs compara- 
ble to typical interstellar absorption might suggest 
that the large Si IV EWs are an artefact of poor 
fitting and/or simply the use of too simplified a 
model. However, for C the situation is different, 
even ignoring the complicated C IV doublet. The 
C II blend is too deep, though just consistent with 
the maximum interstellar EW value, and we also 
identified a number of C I absorption lines in our 
spectrum, of which only C I A1329 is normally 
seen (see e.g. Morton et al. 1972). In the case 
of C IV even trying double-peaked emission had 
very little effect on the amount of absorption re- 
quired, and here the excess is large. The most 
likely explanation would seem to be excess C, at 
various ionization stages, local to the system. In- 
deed, SCH01 found that fits to their X-ray spectra 
were improved by including the effects of a carbon 
K edge with an overabundance of C relative to H 
of (5 — 18) x solar. 

[!htb] 

4-1.2. Emission line strengths 

Once we have taken into account the pro- 
nounced absorption for both the Si IV/O IV] 
blend and C IV with our Gaussian fitting, these 
emission features turn out to be the strongest, 
both with EW^4A. The various O emission lines 
are also generally prominent, the more typical 
lines of O V and O IV having EW~3 and 2 A 
respectively, whilst even the unusual O III] fea- 
ture yields ~1.5A. The weakest line we measure is 
that of C III with only ~0.5A. We note that both 
N VA1240 and He IIA1640, common in a large 
variety of high excitation objects, are very weak 
or entirely absent in our spectrum of 4U 1626-67. 

We also compare our results on the emission 
line strengths of 4U 1626-67 with those published 
for Her X-l (Anderson et al. 1994; Boroson et al. 
1996), Sco X-l (Kallman et al. 1998), and XTE 
J1859+226 and XTE Jl 118+480 (Haswell et al. 
2002) based on other HST data. These systems 
represent respectively, an intermediate mass X-ray 
pulsar, the "prototypical" LMXB, and two rather 
different soft-X-ray transients (black-hole candi- 
dates) in outburst; hence they should provide a 
reasonably diverse context for our study. 

Common emission line characteristics of these 
comparison objects include: (i) very prominent 
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Table 2 

Comparison of 4U 1626-67 Absorption line EWs to Interstellar lines from HST Quasar 

Key Project 



ID 


^vac 


Flux (10~ 16 
ergcm~ 2 s -1 ) 


EW 


EW from HST Key Project 




(A) 


(kms- 1 ) 


Median (kms x ) 


Maximum ( km s 1 ) 


CI(?) 


1188.99 


-3.88 ± 1.61 


67 ±28 






Si II 


1190.42 


-6.53 ± 1.39 


113 ±24 






Si II 


1193.29 


-9.20 ±0.97 


159 ± 17 






S III 


1194.06 


-9.20 ±0.97 


159 ± 17 






NI 


1200.13 


-12.28 ±0.99 


219 ± 18 


240 ± 65 


260 


SII 


1250.58 


-10.49 ±2.80 


182 ± 49 






SII 


1253.81 


-5.91 ±1.32 


101 ±23 


70 ±90 


220 


Si II 


1260.42 


-12.63 ±0.94 


221 ±18 


190 ±65 


320 


j£n> 


1259.52 












1277.41 


-1.89 ±0.79 


31 ± 13 






CI 


1280.33 


-1.21 ±0.99 


20 ± 16 






01 


1302.17 


-11.99 ±2.73 


187 ±43 


105 ± 55 


220 


Si II 


1304.37 


-6.84 ±1.86 


106 ± 29 


95 ±35 


190 


CI 


1328.83 


-2.88 ±1.03 


46 ± 16 






CII 


1334.53 


-16.59 ±3.10 


260 ± 50 


160 ± 45 


230 


+CII 


1335.71 










Si IV 


1393.76 


-10.25 ±1.17 


163 ± 19 


60 ±40 


160 


Si IV 


1402.77 


-13.81 ±1.32 


220 ±21 


40 ± 10 


60 


Si II 


1526.71 


-7.17 ±0.79 


109 ± 12 


95 ±25 


150 


CIV 


1548.20 


-14.21 ±1.91 


220 ± 30 


80 ±25 


110 


CIV 


1550.77 


-9.92 ± 1.79 


153 ± 28 


45 ±20 


80 


Fe II 


1608.45 


-5.87 ± 1.56 


94 ±25 






CI 


1657.22 


-2.51 ±1.09 


42 ±18 






Al II 


1670.79 


-9.48 ±1.52 


164 ± 26 


115 ±30 


220 



a As presented in figure 5 of Savage et al. 2000 
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N V A 1238, 1242 doublet, (ii) equally strong or 
stronger C IV A1548, A1550 doublet (except for 
XTE J1118+480, for which abundance anomalies 
are postulated), (iii) prominent He IIA1641 , (iv) 
similar strength Si IV/O IV] blend, and (v) usu- 
ally slightly weaker O V and C IIIA1175 (again 
except for XTE Jl 118+480). Hence, the absence 
of noticeable N V and He II is one distinctly un- 
usual feature of 4U 1626-67. The lack of He II 
provides strong support for the conclusion of the 
X-ray abundance analysis of SCH01, namely that 
the white dwarf has no He to transfer, but is in 
fact the core of a more massive initial star, com- 
prised largely of either C-O-Ne or O-Ne-Mg. The 
apparent reduction in N abundance is then not 
unexpected. Furthermore, it also appears that 
there is an excess of O in 4U 1626-67; O IVA1342 
is very strong here whereas it is only otherwise 
noted in the case of Sco X-l (although it may 
also be present in the published spectrum of XTE 
J1859+226), whilst O III] multiplet A1661-1669is 
not identified in any of these comparison systems. 

In terms of distinguishing between the two pos- 
sible white dwarf compositions, Mg is the key. As 
SCH01 discuss, either C-O-Ne or O-Ne-Mg white 
dwarf models can account for the high observed 
Ne and C abundances. In the FUV the only Mg 
interstellar absorption line that is sometimes seen 
is Mg II at 1240 A but in the NUV there is the 
Mg II A2796, A2803 doublet, which again appears 
as an interstellar absorption feature but is also 
found in emission in a range of astrophysical situa- 
tions, including quiescent LMXBs (McClintock & 
Remillard 2000). It is possible that there is weak 
A1240 absorption, whilst analysis of the A2796, 
A2803 doublet (in our NUV spectrum), is rather 
complex and inconclusive. The lower spectral res- 
olution of the NUV spectrum makes any disen- 
tangling of emission versus absorption components 
difficult. There is marginal evidence for emission, 
which would in turn imply probable excess ab- 
sorption, but equally if we assume no emission 
the absorption is comparable to what we would 
expect from interstellar. It would seem that fur- 
ther higher signal-to-noise/resolution spectra con- 
centrating on this Mg II doublet would be helpful 
to derive observational constraints from this av- 
enue. 



4.2. Emission region dynamics and loca- 
tion 

In terms of the emission line profiles we find two 
broad categories from our fitting; those consistent 
with a single emission region and those preferring 
a Doppler pair of emission regions. C III and C IV 
can be represented by emission from a single re- 
gion with a broad dispersion in velocities giving 
a FWHM^1500kms _1 for the lines, with a neg- 
ligible shift in line centers. Si IV is similar (con- 
sidering the final fit with double-peaked O IV]), 
as is O III], although the FWHM-gOOkms" 1 and 
600 km s -1 are somewhat smaller and 100 km s -1 
velocity shifts were indicated, though for O III] 
and possibly for Si IV these may well not be 
significant. By contrast, the more highly ion- 
ized O V singlet most likely requires a Doppler 
pair of emission regions (with V ~ —1000 and 
+600 km s -1 ), and similar results arc suggested 
for the O IV triplet. Even then the red-shifted 
components have very large FWHM^2000 kms -1 , 
roughly twice that of the blue-shifted components, 
which may indicate that the there are in fact mul- 
tiple red-shifted emission regions. 

Detailed modeling of the various emission re- 
gions and their dynamics is beyond the scope of 
this paper. However, some physical insights can 
be gained by considering results from earlier gen- 
eral studies of UV line formation. It should be 
noted that all these studies assumed typical cos- 
mic abundances for the emitting gas, in particular 
that the heavier elements have trace abundances 
relative to H and He, which is almost certainly not 
the case for the transferred material in 4U 1626-67. 

Kallman & McCray (1982) presented a range 
of theoretical models for the ionization structure 
of gas irradiated by a powerful X-ray source, cov- 
ering a range of astrophysically relevant param- 
eters. This enables us to look at the differing 
dependence of the various ionic species on the 
ionization parameter (£ = Lx/nR 2 , where Lx 
is the X-ray irradiating luminosity, n is the gas 
density and R is the distance from the emission 
region to the X-ray source). Their 8 models in- 
clude different values of the luminosity, its spec- 
tral shape and the gas density. The case of a 
high density gas (n — 10 n cm~ 3 ), illuminated by 
a 10 keV thermal bremsstrahlung X-ray spectrum 
with Lx = 10 37 ergs -1 with line trapping effects 
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included, is particularly interesting. Such condi- 
tions may well be similar to those in an accretion 
disc corona/atmosphere above the disc of 4U 1626- 
67. The calculations show that the peak of the 
predicted abundances for our two groups of emis- 
sion lines are separated in terms of the value of £ 
required. C III, Si IV and O III all peak at around 
log£ = 1.35, whereas O IV and V require a signif- 
icantly higher log£ ~ 1.65. C IV lies between the 
two in terms of £, but still its abundance is neg- 
ligible by the time you reach the peak log£ value 
for the more highly ionized O species. Although 
the details will not be exactly applicable to our 
case, the important point is that such a division 
in terms of the values of £ is plausible. In con- 
sequence, there could be differences in the radial 
dependence of the emission. For instance, model- 
ing of the broad emission lines of Her X-l (admit- 
tedly for a posited disk wind situation) by Chiang 
(2001) found that the outer radii of the emission 
region for a given ion were systematically smaller 
for the higher ionization lines. Similar "ionization 
stratification" was also found from reverberation 
mapping of the broad line regions of AGNs (see 
e.g. Krolik et al. 1991; Korista et al. 1995). The 
effect is probably quite subtle, and in fact specific 
modeling of the emission lines from X-ray heated 
accretion discs by Ko & Rahman (1994) did not 
find any clear difference in the radial dependence 
of the line formation between C IV and O V, al- 
though for Si IV the contributions did drop more 
rapidly with decreasing radial location. They also 
predicted that overall the greatest contributions to 
the UV lines should be from the outermost radii. 

Observationally, the trend in £ does appear to 
follow that of the velocity spread/FWHM of our 
emission lines. For double-peaked emission lines 
the velocity shift of the peaks should be a measure 
of the Keplerian velocity of the outermost con- 
tributing radial annulus (see Smak 1981; Home & 
Marsh 1986), whilst for apparently single-peaked 
lines the HWHM should also be an approximate 
indicator. We note that both the single Gaussian 
and the twin-Gaussian profiles are merely different 
approximations to the actual velocity profile one 
would expect from gas distributed in a Keplerian 
flow. In fact, the separation into the two groups 
is probably a consequence of which is the better 
approximation for a given velocity spread. 

To aid in the comparison of observations with 



the above models we may also consider the even 
broader X-ray lines. SCH01 found that all these 
lines had similar overall width, and hence could 
be fitted with double-peaked profiles. In their 
case, the FWHM of the blue and red-shifted 
components could be taken as equal at ^2000- 
3000 km s -1 , with large velocity shifts for peaks 
ranging from - 1600 km s -1 to -2600 km s -1 in the 
blue and 800-1900 km s" 1 in the red. As they 
comment, these line profiles could also be con- 
sistent with a disc origin. If one considers the 
velocities of the various X-ray and FUV lines in 
terms of a Keplerian accretion disc or atmosphere 
just above (assuming i = 30°, the upper limit for 
a 0.02M Q donor), one would place the outermost 
contributing radii for C III, Si IV and O III at the 
outer (tidally truncated) edge (r t — 2 x 10 10 cm, 
where Vsini = 530 km s -1 ), C IV a little far- 
ther in, O IV and O V slightly farther in still 
and for the X-ray lines from r ~ 3 x 10 10 cm 
to the inner edge (at magnetospheric corotation, 
r co = 6.5 x 10 8 cm, where Vsini = 2900 kms -1 ). 
But this is where the simple physical picture out- 
lined above appears to break down. As discussed 
by SCH01, given the luminosity 4 of the source 
(L x = 2 x 10 36 (d/3kpc) 2 ergs" 1 ), and a gas den- 
sity of n ~ 10 n cm- 3 requires R> 10 10 (d/3kpc), 
comparable to the radius of the outer disc to ob- 
tain £ < 10 3 erg cm s _1 and allow formation of 
even the more highly ionized X-ray line emitting 
species. For our FUV lines, with £ ~ 50 erg cm 
s , the situation becomes even worse, requiring 
i?>10 n (d/3kpc), almost an order of magnitude 
larger than the entire 41.4 min period binary! 
One solution might be to assume a higher density; 
indeed densities exceeding 10 16 cm -3 may be pos- 
sible in the accretion disc proper. In any case, we 
suggest that the velocities are more reliably indi- 
cating the radial limits of the emission locations, 
and that it is the details of the ionization structure 
that need further scrutiny. Admittedly, there are 
outstanding issues related to the velocity struc- 
ture too, notably the asymmetry of the widths of 
the Doppler pairs for O V (or an extra redshifted 
component), and for all the pairs of both X-ray 
and FUV lines the inequality of the red and blue 
velocity shifts (Vbiuc — 2 x U rc d)- 

4 A minimum distance of 3 kpc was derived by Chakrabarty 
(1998) based upon the mass transfer rate implied by the 
spin-up rate of the pulsar in 4U 1626-67. 
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4.3. Conclusions 

Analysis of a 32ks timc-averagcd FUV spec- 
trum of the ultra-compact X-ray binary, 4U 1626- 
67 has provided further observational constraints 
on the nature of the donor star and the kinemat- 
ics of the line emission regions in the system. We 
find evidence for both a lack of He and N, to- 
gether with an excess of O in the line emitting 
gas. Moreover, additional C absorption local to 
the system is probable. These observations are 
fully consistent with the X-ray spectroscopic re- 
sults of Schulz et al. (2001), which indicated that 
the donor is a very low-mass C-O-Ne or O-Ne- 
Mg white dwarf. In addition, we also find com- 
plicated velocity structures in emission, similar to 
the Schulz et al. Doppler pairs of X-ray lines. All 
the FUV lines are very broad, with the highest 
ionization parameter species exhibiting possibly 
double-peaked structures. The interpretation of 
these line profiles as a consequence of emitting re- 
gions in a Kcplcrian disc about the neutron star 
is certainly plausible, though the simplest models 
may not be adequate. In any case, the combina- 
tion of high quality FUV and X-ray spectral line 
data may provide modelers with a rich data set 
for follow-on detailed studies of the disk dynamics 
and ionization structure. 

Support for this work was provided by NASA 
through STScI grant GO-06624.01 and LTSA 
grant NAG5-7932. 
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